High-field magnetization study of the S = 1/2 antiferromagnetic Heisenberg chain 
[PM Cu(N0 3 )2(H 2 0)2]„ with a field-induced gap 
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We present a high-field magnetization study of the 5=1/2 antiferromagnetic Heisenberg chain 
[PM Cu(N03)2(H20)2]n- For this material, as result of the Dzyaloshinskii-Moriya interaction and a 
staggered g-tensor, the ground state is characterized by an anisotropic field-induced spin excitation 
gap and a staggered magnetization. Our data reveal the qualitatively different behavior in the di- 
rections of maximum and zero spin excitation gap. The data are analyzed via exact diagonalization 
of a linear spin chain with up to 20 sites and on basis of the Bethe ansatz equations, respectively. 
For both directions we find very good agreement between experimental data and theoretical calcu- 
lations. We extract the magnetic coupling strength J jks along the chain direction to 36.3(5) K and 
determine the field dependence of the staggered magnetization component m s . 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.30.Gw, 75.50.Xx 



Motivated by the rich variety of different magnetic 
ground states, such as quantum critical behavior or gaps 
in the spin excitation spectra, quasi-one-dimensional 
quantum magnets have been the focus of intense exper- 
imental and theoretical research efforts in recent years 
iSiiii. To gain deeper insight into the physics of such 
quantum spin systems well-defined model compounds 
need to be explored. Here, the uniform 5=1/2 anti- 
ferromagnetic Heisenberg chain (AFHC) is of particular 
interest, since it is exactly solvable using the so-called 
Bethe ansatz equations 5 i 6 i 7 . 

In S — 1/2 AFHCs, lacking inversion symmetry, ad- 
ditional terms in the Hamiltonian have to be taken into 
account, that is the Dzyaloshinskii-Moriya (DM) interac- 
tion and an alternating g tensor This gives rise to an 
effective staggered field h s perpendicular to the applied 
magnetic field H. Then the Hamiltonian is written as^i 

H = J^[S,S l+1 - h u Sf - (-iyh s Sf] (1) 

i 

with J as the coupling constant, h u — g^sH/J as the 
effective uniform field, and h s the induced effective stag- 
gered field. In the following we refer to this as the stag- 
gered S — 1/2 AFHC model. Resulting from this exten- 
sion of the uniform 5=1/2 AFHC are the opening of an 
anisotropic spin excitation gap with application of a mag- 
netic field and new, particlelike excitations such as soli- 
tons, antisolitons, and their bound state, the "breather" 
22*il. Moreover, by fully evaluating the effect of the DM 
interaction on the ground state properties, a crossover 
to a qualitatively different high-field behavior has been 
predicted recently—. 



The model for the staggered 5=1/2 AFHC has 
been used to describe two materials in particular, 
copper benzoat© 3 " and copper pyrimidine nitrate [PM 
Cu(N03) 2 (H 2 0)2]n — • For the latter compound, from 
a single-crystal study a magnetic exchange parameter 
J/ks = 36 K is derived. Further, an additional Curie-like 
contribution to the magnetic susceptibility at low tem- 
peratures is observed, which varies strongly with magni- 
tude and direction of the applied external field. Specific- 
heat measurements in magnetic fields verify the predicted 
formation of an anisotropic spin excitation gap, whose 
magnitude also depends on size and orientation of the ex- 
ternal field—. The spin excitation gap and the Curie-like 
contribution to \ are largest for the same field direction 
(referred to as c" , for notation see Ref. 13) and vanish 
for one direction perpendicular to c" in the a — c plane 
(referred to as a"). 

For 5 = 1/2 AFHC materials, (high-field) magneti- 
zation experiments are abundant and are perfectly de- 
scribed by theory ii^iSiiSiiLiS. Recently, for the stag- 
gered 5 = 1/2 AFHC the magnetization curve has been 
calculated by several groups iiisLSfliSI. As yet, these 
theoretical predictions have not been verified experimen- 
tally. 

Therefore, in this Rapid Communication we present a 
magnetization study on [PM Cu(N03)2(H 2 0)2] n cover- 
ing the entire field range up to saturation, i.e. fioH = 53 
T. With our study of the magnetization along a" and c" 
we establish the contrasting behavior along these two di- 
rections, the first representing the uniform, and the latter 
the staggered 5=1/2 AFHC. The behavior of the uni- 
form 5=1/2 AFHC is evaluated on basis of the Bethe 
ansatz equations iA£. In contrast, for the staggered 5 



= 1/2 AFHC we analyze our data by means of exact di- 
agonalization of linear chains with up to N = 20 spins, 
based upon the staggered field theory by Oshikawa and 
Affleck 8 - 9 . From this analysis we find very good agree- 
ment between experimental data and theoretical calcula- 
tions. We determine the characteristic parameters, i.e., 
the coupling constant J/ks and the staggered magneti- 
zation m s . 

Comparing our finite-size calculations to previous 
density-matrix renormalization group (DMRG) studies 
^2ii2i22i2i for the Hamiltonian Q , we find perfect agree- 
ment between our results and those of other groups. 
Only for low- fields finite-size effects are present. The 
advantage of our finite-size calculations compared to the 
DMRG method lies in its simplicity and short computa- 
tion time. 

Single crystals of [PM Cu(N0 3 ) 2 (H 2 0) 2 ]„ have been 
grown by slow evaporation of the equimolar aqueous so- 
lution of copper nitrate and pyrimidine— . The crystals 
show well-defined facets and the principal axes can be 
identified easily. We have checked by low-field magne- 
tization measurements that the magnetic susceptibility 
matches the one published in Ref. 13. For the magne- 
tization measurements the samples were oriented along 
the characteristic orientations a" and c" (misalignment 
< 5°), glued to the tip of a plexi glass rod and placed 
inside a thin walled teflon cylinder. The magnetization 
signal of the sample holder was negligible. 

Magnetization measurements were carried out at the 
Laboratoire National des Champs Magnetiques Pulses in 
Toulouse in pulsed magnetic fields up to /j,qH = 53 T. 
Pulsed magnetic fields were obtained by discharging a 
capacitor bank in a solenoid according to a crowbar de- 
scribed in Ref. 23. The pulse duration was about 200 
ms with an increasing time of 25 ms. The magnetization 
was detected as a voltage V induced in a compensated ar- 
rangement of pick-up coils wound concentrically around 
the sample and coupled to it with the coupling constant 
rj, such that V = rjflSM/St, with as the sample volume 
— . The absolute magnetization was obtained by numeri- 
cal integration of this voltage. Due to the limited sample 
space (diameter < 1.6 mm), the absolute signal was small 
(< 10~ 5 Am 2 ). To achieve a higher accuracy of the sig- 
nal calibration additional measurements were performed 
in magnetic fields up to 5 T in a commercial supercon- 
ducting quantum interference device magnetometer. 

In Fig. 1 we present the magnetization curve of [PM 
Cu(N0 3 ) 2 (H 2 0) 2 ]„ as a function of field at 1.6 and 4.2 
K for the two characteristic orientations, i.e., H\\a" (Fig. 
1(a)) and \\c" (Fig. 1(b)), respectively. Comparing the 
magnetization along the two directions, an anisotropic 
response is observed. At 1.6 K for H\\a" we find the 
archetypical behavior of the 5=1/2 AFHC — . In con- 
trast, for fields parallel to the c" axis an additional low 
field contribution and a delayed saturation of the mag- 
netization occurs. To emphasize this difference in the 
insets of Figs. 1(a) and (b) we plot the derivatives of 
the magnetization, dM / d(fioH) . For small fields ||c" the 
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FIG. 1: The field dependence of the magnetization of [PM 
Cu(NC>3)2(H20)2]n with the external field aligned along the 
a" (a) and the c direction (b). In the insets, the field deriva- 
tive dM/d(/j,oH) is displayed. 

initial slope is more than twice as large as for the a" axis. 
At high fields (> 35 T), the saturation of the magneti- 
zation for i?||c" is suppressed, as indicated by a smaller, 
broader feature in dM/d(fj,oH), as compared to the a" 
axis response. Increasing the temperature to 4.2 K re- 
duces the difference in the M{hqH) curves between the 
two directions, but does not completely suppress it. 

The deviation from the uniform 5=1/2 AFHC behav- 
ior along the c" direction is attributed to an additional 
magnetization component. It increases much faster and 
passes through a maximum at a much lower field than 
the uniform saturation field noH sat . Since the second 
component is not present along a" , we ascribe it to the 
staggered magnetization m s . 

The g tensor of [PM Cu(N0 3 ) 2 (H 2 0) 2 ]„ has been de- 
rived from electron spin resonance measurements—. In 
the uniform 5=1/2 AFHC model the saturation field is 
calculated according to the formula^ 

H c = AJS/g^B- (2) 

For zero temperature and using J/fcs =36±0.5K, 
g a " = 2. 14 ±0.02, and g c » = 2.21 ±0.02 from Ref. 13, we 
obtain fi H sat = 50.1 ±0.8T and 48.5 ± 0.8 T along a" 
and c" axes, respectively. The saturation magnetization 
M s is calculated to M Sja » = 1.07 ± 0.01 /ib/Cu atom and 



3 



M s>c " = 1.11 ± 0.01 hb/Cu atom. Thus, for the uniform 
5 = 1/2 AFHC at T = 1.6 K < J/k B , the saturation 
magnetization should be approached at highest experi- 
mental applied field. Indeed for H\\a" (Fig. 1(a)), the 
T = 1.6 K curve has an initial slope lower than at 4.2 
K. With increasing field the curvature becomes larger, 
crosses the 4.2 K curve near 38 T and almost reaches 
saturation at ss 53 T. The data for both temperatures 
become nonlinear with field for fioH > 15 T. Moreover, 
with decreasing temperature the data sets approach the 
T = curve for the uniform 5 — 1/2 AFHC in full agree- 
ment with previous experimental work—. 

From a theoretical point of view, the magnetization 
curve of the uniform 5=1/2 AFHC has been computed 
for T — by Bethe ansatz — . However, our measure- 
ments have been carried out at temperatures T > 0.04 J, 
where thermal fluctuations lead to a substantial round- 
ing of the cusp at the saturation field in the T — mag- 
netization curve. On the other hand, at 1.6 K we have 
T < 0.05J, substantially smaller than the lowest cor- 
responding temperature studied in Ref. 17 for copper 
pyrazine dinitrate. Even at T = 0.05J, the magnetiza- 
tion curve computed for the uniform 5=1/2 AFHC on a 
ring with N = 16 sites still exhibits clear finite-size ef- 
fects, precluding an analysis along the lines of Ref. 17. 
We therefore use results obtained by the thermodynamic 
Bethe ansatz in the thermodynamic limit N = oo and at 
arbitrary temperature T to describe the magnetization 
curve for i? || a"i£. 

The situation for H \\ c" is quite different. On the one 
hand, the Hamiltonian (JIJ cannot be solved exactly with 
a nonzero staggered field h s and we therefore have to rely 
on a numerical treatment. Whereas the ^-component 5 2 
of the total spin is conserved for h s = 0, even this is 
not the case anymore for h s ^ 0. This has two conse- 
quences: («) The reduced symmetry of the Hamiltonian 
restricts the system sizes N that can be accessed, and 
(ii) each new set of magnetic fields h u and h s requires a 
new numerical determination of the ground state. On 
the other hand, the field-induced opening of a gap 
leads to the following two simplifications: (i) For the 
high magnetic fields studied here, the gap is sufficiently 
large to suppress thermal excitation at low temperatures. 
Therefore, finite temperature is expected to have only a 
comparatively small effect. Indeed, this is confirmed by 
the difference of the T — 1.6 K and 4.2 K curves in Fig. 1, 
which is noticeably smaller along the c" direction (panel 
(b)) than for the a" one (panel (a)). This permits us to 
compare a measurement at low but finite T with a com- 
putation at T = 0. (ii) In most field ranges considered 
here, the correlation length turns out to be sufficiently 
short such that finite-size effects can be neglected already 
for systems with only N = 20 sites. More precisely, the 
correlation length is large only for a small staggered field 
h s and thus only the low-field region suffers from finite- 
size effects. These lead to an artificial low-field peak in 
the staggered magnetization m s (see e.g., dashed curve 
in Fig. 2(b)) whose position roughly determines the re- 
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FIG. 2: Experimental and theoretical magnetization curves 
for [PM Cu(N0 3 )2(H 2 0)2]n at T = 1.6 K for different field 
directions, (a) HoH\\a": solid line - experimental data, dashed 
line - fit assuming uniform S = 1/2 AFHC for T = 1.6 K, dot- 
ted line - T = calculations—-, (b) fj,oH\\c": dotted line - 
uniform magnetization m u , dashed line - staggered magneti- 
zation m s , dash-dotted line - calculated physical magnetiza- 
tion niphys for c = 0.11, solid line - experimental data, (c) 
fM)H\\c": dashed line - m p hy S for c = 0.11, dotted line - m p h V s 
for c = 0.2, solid line - experimental data. 

gion up to which finite-size effects are still relevant, as 
evidenced by comparison with results for N < 16 (not 
shown). Due to the fast disappearance of finite-size ef- 
fects for higher magnetic fields it is completely sufficient 
for our purposes to apply the Lanczos diagonalization 
procedure to rings with N < 20 sites. Therefore, the 
additional effort of a DMRG procedure 12 i 19 i 20 i 21 i s no t 
necessary here. 

A final remark is in order before we present our nu- 
merical results and compare them with the experiment. 
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In the presence of a staggered field h s = ch u , which is 
related by a constant anisotropy parameter c to the uni- 
form field h u , the physical magnetization m p h ys (which 
is measured experimentally) is given by the superposi- 
tion of the uniform m u and staggered m s magnetization 
components 

m p hys = m u + cm s . (3) 

To our knowledge, previous numerical works 
only showed m u and m s separately, while numerical re- 
sults for the combination m p h ys (Eq. EJ have not been 
published yet. 

For the magnetization along a" (Fig. 2(a)), using 
fioHsat — 50.6 T, g a n = 2.14, and taking into account the 
finite experimental temperature, we find very good agree- 
ment between our data and the Bethe ansatz result — . 
The deviations for fields ^i n H > 50 T can be attributed 
to the misalignment of the crystal, while for smaller fields 
experimental data and calculated result match within 
2%. Thus, for this field direction the uniform S = 1/2 
AFHC is established. 

In Fig. 2(b) we depict the field dependence of the 
magnetization along c". For fields noH > 10 T we can 
fully describe these data on basis of the staggered 5 — 
1/2 AFHC model using ^ H sat = 49.3 T, g c „ = 2.19. 
Here, we have used an anisotropy parameter c = 0.11. 
This value has been obtained from a comparison of the 
data and calculations for c £ [0.08;0.12], step size 0.01, 
as the optimum solution. The values /j,qH along a" and 
c" are fully consistent with those predicted from Eq. (2) 
within their error bars. On average, the saturation fields 
1 1 c" and 1 1 a" correspond to a magnetic coupling strength 
J/k B = 36.3(5) K. 

With the value of the anisotropy parameter c we can 
decompose m p h vs into the uniform (ra u ) and staggered 
(to s ) component, according to Eq. 10. Both, m u and 
m s , are included in Fig. 2(b). Their field dependence 
closely resembles the ones obtained in Ref. 20 for the 
case h u — 10 h s . Specifically, we find that m s traverses 



a maximum at ~ 40 T, while m u and m s approach finite 
but nonsaturated values for largest fields. Our analysis 
establishes the staggered 5=1/2 AFHC for the c" axis 
of [PM Cu(N0 3 ) 2 (H 2 0) 2 ]„. 

Further, we performed calculations for c values in the 
range 0.08-0.28. The cases c= 0.2 and 0.11 are depicted 
in Fig. 2(c). From these data it appears that for decreas- 
ing parameter c the curvature of m p hy S increases and 
approaches the uniform S — 1/2 AFHC magnetization 
behavior for vanishing staggered field h s , as expected. 
With the definition of the Hamiltonian of the staggered 
5=1/2 AFHC model in Ref. 13 the anisotropy parame- 
ter c = 0.11 from this work translates into a value 0.24 25 . 
Thus, we find perfect agreement of our c value with the 
one obtained from magnetic susceptibility measurements 
in Ref. 13, i.e., c = 0.235. 

In conclusion, we have performed high-field magnetiza- 
tion experiments on [PM-Cu(N03)2-(H 2 0)2]n at temper- 
atures 1.6 and 4.2 K. We analyzed our data by means of 
Bethe ansatz equations for the experiments along the a" 
axis and exact diagonalization of linear chains containing 
N = 20 spins, based upon the staggered field theory by 
Oshikawa and Affleck for the c" axis. The very good 
agreement of our data with our theoretical calculations 
and those of other groups iiL2& verifies the predictions 
for the uniform and staggered 5=1/2 AFHC models. 
This way, for the staggered case we have extracted the 
staggered magnetization component m s (fioH). 

Recently, for the staggered 5=1/2 AFHC model a 
low-field/high-field crossover in the staggered magnetiza- 
tion has been predicted^. For [PM-Cu(N0 3 )2-(H 2 0)2]„ 
these calculations would imply deviations from the stag- 
gered behavior in the field range fioH > 50 T, and thus 
cannot be observed in our study. To verify the predic- 
tions from Ref. 12, in the future it will be interesting 
to perform analogous magnetization studies on related 
materials with smaller J values. 

This work has partially been supported by funds of the 
European contract no. HPRI-CT-1999-40013 and by the 
DFG under contract no. SU229/6-1. 
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